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We investigate the quark orbital angular momentum of the nucleon in the absence of gauge- 
field degrees of freedom, by using the concept of the Wigner distribution and the light-cone wave 
functions of the Fock state expansion of the nucleon. The quark orbital angular momentum is 
obtained from the phase-space average of the orbital angular momentum operator weighted with 
the Wigner distribution of unpolarized quarks in a longitudinally polarized nucleon. We also 
derive the light-cone wave function representation of the orbital angular momentum. In particular, 
we perform an expansion in the nucleon Fock state space and decompose the orbital angular 
momentum into the iV-parton state contributions. Explicit expressions are presented in terms of 
the light-cone wave functions of the three-quark Fock state. Numerical results for the up and down 
quark orbital angular momenta of the proton are shown in the light-cone constituent quark model 
and the light-cone chiral quark-soliton model. 

PACS numbers: 12.38.-t,12.39.-x,14.20.Dh 



1 



I. INTRODUCTION 



The spin structure of the nucleon is one of the most important open questions in hadronic 
physics which has recently attracted an increasing attention. While the quark spin contri- 
bution is well known from the accurate measurements in polarized deep inelastic scatter- 
ing l|, |2J, the situation is much more unclear for the remaining contributions from the quark 
orbital angular momentum (OAM) and those from the gluon spin and OAM. In the present 
paper we will focus on the quark OAM. As a first step, we ignore the contributions from the 
gauge-field degrees of freedom. Under this assumption, there is no ambiguity in the defini- 
tion of the quark OAM operator 3r|5| . To unveil the underlying physics associated with the 
quark OAM, we discuss two different representations. The first one is obtained by using the 
concept of the Wigner distributions jfj], whereas the second one is based on light-cone wave 
functions (LCWFs). The extension to include the gauge- field degrees of freedom is left for 
a further study. 

The Wigner distribution was originally constructed as the quantum mechanical analogue 
of the classical density operator in the phase space. Recently, it was also exploited to provide 
a mult i- dimensional mapping of the partons in the nucleon j§-8|. An important aspect of 
the Wigner distribution is the possibility to calculate the expectation value of any physical 
observable from its phase-space average with the Wigner distribution as weighting factor. In 
particular, we will show that the quark OAM can be obtained from the Wigner distribution 
for unpolarized quarks in the longitudinally polarized nucleon. 

On the other hand, the LCWFs provide the natural framework for the representation of 
OAM. This is due to the fact that the LCWFs are eigenstates of the total OAM for each 
iV-parton configuration in the nucleon Fock space. In particular, we will consider the three 
valence-quark sector, giving the explicit representation of the quark OAM in terms of the 
partial-wave amplitudes of the nucleon state. 

The plan of the paper is as follows. In Sec. Hi] we give the relevant definitions for the 
quark OAM operator and derive the expression of the OAM in terms of a Wigner distri- 
bution. In Sec. II III we first discuss the LCWF representation of the quark OAM, giving 
the decomposition into the contributions from the 7V-parton states. For the three-valence 
quark sector we further derive in Sec. |IV] the explicit expressions for the contributions from 
the different partial-wave amplitudes of the nucleon state to the OAM. The corresponding 
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expressions for the distribution in x of the OAM are collected in the Appendix. The for- 
malism is finally applied in a light-cone constituent quark model and the light-cone version 
of the chiral quark-soliton model, and the corresponding results are discussed in Sec. [V] 
Concluding remarks are given in the last section. 



II. QUARK ORBITAL ANGULAR MOMENTUM 
A. Definitions 

If one neglects gauge-field degrees of freedom, the quark OAM operator for a given flavor 
q can unambiguously be defined as 

L\ = J dr~ d 2 r^p 9 (r-,r)-f + (rx V r ) ^(r",r), (1) 

where normal ordering is understood and the transverse vector a = (a 1 , a 2 ) has been in- 
troduced for the generic 4-vector a. Furthermore, in Eq. (JT]) we used for convenience the 
notation V=V — V. Using the Fourier transform of the quark field from the coordinate 
space to the momentum space 

^e-^-H^^^), (2) 



the quark OAM operator can alternatively be written as 
We define the quark OAM l\ as the following matrix element of the quark OAM operator 



/ ^?^( fc+ » fe )7 + (fe x V fe )/ 9 (fc + ,fe). (3) 



(p',+\Ll\p,+)=£l(p',+\p,+), (4) 
where the momenta of the incoming and outgoing nucleon are given by 

P = P-#, p' = P + f- (5) 
Since the nucleon states \p, A) in Eq. (j4j) are normalized as 

+\ Pi +) = 2P+ 5 (A+) (2tt) 3 6® (A) , (6) 

the forward limit p' = p has to be treated with care. One can easily get rid of this normal- 
ization by integrating over A + and A 

avoiding in this way to use normalizable wave packets or infinite normalization factors. 



B. Orbital angular momentum from Wigner distributions 



In this section we derive the link between the quark OAM and the Wigner distributions. 
In particular, we show how it is possible to express Eq. (j7j) as the phase-space average of 
the intuitive OAM operator (r x k) weighted with the Wigner distribution of unpolarized 
quarks in a longitudinally polarized nucleon. 

From the definition of the OAM introduced in the previous section, it is natural to 
interpret the integrands of Eqs. ([T|) and ([3]) as the quark OAM density operators in position 
and momentum space, respectively, 

L q z (r~, r) = — ip 9 (r~, r)^ + (rx V r ) tp q (r~, r), (8a) 

2 \ / 2 

Ll(k + ,k) = -^\k+,k) 1 + (hV^'(fc + ,fc). (8b) 

The generalization to the joint position and momentum space corresponds to the quark 
OAM density operator in the phase space 

Ll(r-,r,k + ,k) = - t - J ^e'^-'^^-r^lrxV^'lr:,^, (9) 

where r, = r + z/2 and rj = r — z/2. This operator contains the OAM density operators in 
position and momentum space as specific limits, i. e. 

J dk + d 2 kL q z (r~,r,k + ,k) =L q z (r',r), (10a) 
J dr~d 2 rL q z (r-,r,k + ,k) =Ll(k + ,k). (10b) 



Introducing the Wigner operator 



1 /"rl7"fl 2 2 

W^ q (r~,r,k\k) = -j ^^e i ( k+z -- k -*H\rJ,r f ) 1 +r(r7,r i ), (11) 
the OAM density operator in the phase space can also be expressed as 

L q z (r~, r, fc+ k) = 2 (r x fc) W [j+]q (r~, r, k + , k). (12) 
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Using the definition of the OAM in Eq. (jTJ), we can write 
dA+ d 2 A 



2P+ 


(2tt) 3 


dA+ 


d 2 A 


P+ 


(2tt) 3 


dA+ 


d 2 A 



(p',+\ J dr~d 2 rdA; + d 2 k L q z {r', r, fc+ k)\p, +) 
dr~ d 2 r dk + d 2 k (r x k) z (p, + \W h+]q {r~ , r, k + , k)\p, +) 
d 2 r dk + d 2 k (r x k) z (p', + \W h+]q (0, r, k + , k) \p, +) J dr~ e iA+r 



P+ (2tt) 3 

= J dxd 2 kd 2 r(r x k) z p [ jl ]q (r,k,x), (13) 

where in the last line we introduced the Wigner distribution for unpolarized quark in a 
longitudinally polarized nucleon p+ + ' 9 ' 

d 2 A 



(P+ f , A'\W^ q (0, r, xP + , k)\P + , -f , A) 



(14) 



with x = k + /P + the fraction of quark longitudinal momentum. Since we consider a nucleon 
state with its transverse center of momentum at the o rigin of the axes, we may identify the 



transverse coordinate r with the impact parameter b 

In the same way as the impact-parameter dependent parton distributions are two- 
dimensional Fourier transforms of the generalized parton distributions (GPDs), the Wigner 
distributions are two-dimensional Fourier transforms of the so-called generalized transverse- 



momentum dependent parton distributions (GTMDs 

d 2 A 



ay, 



13| 



p [ JJl(b,k,x) 



*- b W [ ^(A,k lX ) 



with the GTMD correlator given by 

WS(A, fc, x) = (P + , f , A'\W^ q (0, 0, xP + , k) |P + , -f , A) 



(15) 



(16) 



for the generic twist- 2 Dirac matrix T = 7 + , 7 + 7s, ia^ +r y^ with j = 1,2. The hermiticity 
property of the GTMD correlator 



wff(A,fc,x) 



Wll q (-A,k,x) 



(17) 



implies that p++ 9 is a real quantity, in accordance with its interpretation as a phase-space 
distribution. 

Collecting the main results of this section, we can express the OAM in terms of Wigner 
distributions as 

(18) 



£ q z = J dxd 2 kd 2 b(b x k) z p [ J + + ]q (b,k,x), 



or, equivalently, in terms of GTMDs as [6] 
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(1 



dx d 2 k 



i(kxV A ) z W^ ]q (A,k,x) 

r k 2 

- J dxd 2 k — Ff 4 (x,fc 2 ,0,0), 



A=0 



where F± 4 follows the notation of Ref. 13]. 



(19) 
(20) 



III. LCWF OVERLAP REPRESENTATION 

Following the lines of Refs. HQ , we can obtain an overlap representation of Eq. ( j20j) 
in terms of LCWFs. Since the quark OAM operator is diagonal in light-cone helicity, flavor 
and color spaces, we can write the quark OAM as the sum of the contributions from the 
iV-parton Fock states 

£2 = £^, (21) 



N,/3 



where 



N 



[dx] N [d 2 k] N J2S m [(fc, x V A ) 2 ^(^^(r* , 



(22) 



i=l 



Note that this expression is consistent with the wave-packet approach of Ref. [if]], except for 
the difference by a factor of 1/2 which follows from different conventions for the light-cone 
components. In Eq. (1221) . the integration measures are given by 



[dx} N 
[d 2 k] 



N 



N 



~[dx 

d 2 k: 



,1=1 

' N 

n 2(2^ 



N 
i=l 



(23) 
(24) 



,i=i 



and the iV-parton LCWF ^^g{r) depends on the hadron light-cone helicity A, on the Fock- 
state composition (3 denoting collectively quark light-cone helicities, flavors and colors, and 
on the relative quark momenta collectively denoted by r. In particular, for the active parton 
i one has 



rr m — rp . 



K 

k out 



ki-{\ 
ki + (l 



X i 



7 



A 

2 ' 

A 

2 ' 



(25) 
(26) 
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and for the the spectator partons j ^ i, one has 



X4 — X 7 - kj — kj -\- Xj 2 ; (^^) 



j ~ 3i 



out 

X j — J J 



k° ut — kj xj 2 . (28) 



Acting with Va on the LCWFs in Eq. ( 1221) and then setting A = leads to the following 
expression 1 

/N N 
[d 2 k] N J2 S mJ2^- x n) [*3j(r) (fc,x V fen ) ^(r)] . (29) 
i=l n=l 2 

This expression coincides with our intuitive picture of OAM. Indeed, since the transverse 
position r is represented in transverse- momentum space by iVfc, it means that 

N N 
n=l n=l 

represents the transverse position of the active quark relative to the transverse center of 
momentum R = ^ n x n r n . Moreover, by construction, fcj represents the relative transverse 
momentum of the active quark. It follows that Eq. ( 129]) gives the intrinsic quark OAM 
defined with respect to the transverse center of momentum 2 . 

An important property of the LCWFs is that they are eigenstates of the total OAM 

N 

- * £ (fc„ x V kn ) z *y r) = l z v^(r) (31) 

71=1 

with eigenvalue l z = (A — A n )/2. As a consequence, from the overlap representation ( 129]) 
it is straightforward to show that the total OAM is given by 



52l*Pm„ (32) 



N. 



1 We find the same expression starting from the representation of the quark OAM operator in momentum 
space given by Eq. ^ and integrating over A to kill the S function subjected to 

2 This is completely analogous to non-relativistic mechanics where the intrinsic OAM is defined relative to 
the center of mass 

^2 X Pn = ^2(r n ~R)X {p„ - X n P) +B.X ^(fn ~ X n P) + ^ X n r n X P 
n n n n 

i 

with R = ^2 n x n f n , P = ^2 n p n , Xn = m n /M and M = X) ra TO ™- This analogy comes from the galilean 
symmetry of transverse space on the light cone. 
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where 



Pm. = £*w; / IMn [#k] s \Kp(r)\ 2 (33) 



is the probability to find the nucleon with light-cone helicity A = + in an iV-parton state 
with eigenvalue l z of the total OAM. 

IV. PARTIAL-WAVE DECOMPOSITION OF QUARK ORBITAL ANGULAR 
MOMENTUM 

In this section we give the explicit expressions for the OAM i\ in terms of the different 
partial-wave components of the three-quark state of the nucleon. 



Working in the so-called "uds" basis 17|, |l8[, the proton state is given in terms of a 



completely symmetrized wave function of the form 

|P, +) = |P, +) uud + |P, +) udu + |P, +) duu . (34) 

In this symmetrization, the state |P, +) u du is obtained from \P,+) UU d by interchanging the 
second and third spin and space coordinates as well as the indicated quark type, with a 
similar interchange of the first and third coordinates for \P,+) duu . For the calculation of 
matrix elements of one-body operators, as in the case of the OAM operator, it is sufficient 
to use only the uud order. As outlined in the previous section, the LCWF ^x'^Tl °f ^ e 
uud component is eigenstate of the total OAM, with eigenvalues l z — (A — J^ i=1 Aj)/2. In 
particular, for a nucleon with helicity A = + we can have four partial waves with l z = 0, ±1, 2 
corresponding to combinations of ^lAaAs w hh appropriate values of the quark helicities Aj. 
However, the 16 helicity configurations of ^^™^ 3 are not all independent. Parity and isospin 
symmetries leave only 6 independent functions ip^ of quark momenta. In particular, the 



complete three-quark light-cone Fock expansion has the following structure 



191421] 



P, +) = |P, +)'«=0 + |p, +)««=! + |p, +)«.="! + |p, + )»»=2 j (35) 
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where 



x ^t4f(i) («}-(2)4+(3) - 4-( 2 >i+( 3 )) ( 36a ) 

| P +) (.=i = /" ^Egs [^^(3)(l,2,3) + ^^ (4) (l,2,3)] 

x ^ (t4f (l)«}-(2)4_(3) - 4f (1)«}-(2)«L(3)) |0>, (36b) 
1^+)'-=-!= f ^EI (-fc 2L )^ (5) (l ; 2,3) 

7 y/XlX 2 X 3 

x (^+(2)4 + (3) - d] + ( 2 )4 + (3)) |0>, (36c) 



7 V /Xi2;2^3 



x (4-(2)4-(3) " «}-(2)4-(3)) |0>- (36d) 

In Eqs. (I36ap -( l36d|) a, (3 = 1,2 are transverse indices, u\ x (ui\) and d\ x (rf iA ) are creation 
(annihilation) operators of u and d quarks with helicity A and color i, and ijjv) are functions 
of quark momenta, with argument i representing Xi and ki and with a dependence on the 
transverse momenta of the form k^ ■ kj only. We also used the notation k^i = k x ± %k y . 
The relations between the wave function amplitudes ^ an d the three-quark LCWFs 



^+^ = ^(i.2)( 1>2 ,3), (37a) 

*^ = V (1,2) (2,1,3), (37b) 

^+|««d = _ ^(1,2) ^ 3 2) + ^(1,2) ^ 3) x) j 5 (37c) 

*+^ = ^ (3,4) (l,2,3), (37d) 

^+^ d = ^ 3 > 4 )(2,l,3), (37e) 

^+|u«d = _ ^(3,4)( 1)3) 2) + ^(3,4)(2 ) 3 ; i)] , ( 37 f) 

tfjjf = -* aL ^ 5 )(l,2,3), (37g) 

*l!!^ = -A; 1Jl A; 3 flV (6) (l,2,3), (37h) 
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where we defined 



V> (1 ' 2) (1, 2, 3) = 2, 3) + ze a % Q fc 2/3 ^(1, 2, 3)] 

^(1,2,3) = [W 3) (l,2,3) + ^ (4) (1,2,3)] . 



(38a) 
(38b) 



The corresponding expressions for hadron helicity A 
parity symmetry 



are obtained thanks to light-cone 



^ -\' 1 -\ 2 -x 3 ({ x ^ kix, hy}) — \& \' 1 \ 2 \ 3 ({xi, ki X ,ki y }). 



A:uud 



(39) 



Using the partial- wave decomposition of the nucleon state in Eqs. (I36al) - fl36dl) we can 
seperately calculate the results for £ q z ' lz corresponding to the contribution of the the quark 
with flavor q in the Fock state component with OAM l z . Summing over all flavors we find, 
in agreement with Eq. (1321) , i l * = l z pi z (we omit the index N = 3 since we considered only 
the three-quark Fock contribution) with p iz = lz (P, +\P, +) lz . In the Appendix we also give 
the results for the partial-wave decomposition of the distribution in x of the OAM for the 
u- and (i-quark contributions. 

For the l z — component, we find: 
• for total w-quark contribution 

C'* = ° = J H 3 [d 2 fc] 3 1 [(1 - xi)k! ■ k 2 + x 2 k\ 

- 2, 3)^ (2) (1, 2, 3) + ^ (1) (2, 1, 3)^(2, 1, 3)" 

^(1,2, 3) + ^(3,2,1) 

((1 - Xl )k! ■ k 2 + x 2 kA^ 2 \l, 2, 3) + (x 2 k x ■ k 3 - z 3 fci • fc 2 )^ {2) (3, 2, 1) 
+ (fei x k 2 ) [^(l, 2, 3) (k x x Vi + k 2 x V 2 )^ (2) (l, 2, 3) 

- V (2) (l, 2, 3) (k! x Vi + k 2 x V 2 V 1} (1, 2, 3) 

- 3, 2) (A* x Vi + fc 2 x V a ) (V 2) (l, 3, 2) - ^ (2) (2, 3, 1)) 
+ ^ (2) (1, 3, 2) (fci x Vi + fc 2 x V 2 ) 3, 2) + ^ (1) (2, 3, 1))] } ; (40) 
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for total (i-quark contribution 
- [dx] 3 [d 2 fc] 



^(1,2,3)^(1,2,3) 



(41) 



xik 2 ■ k 3 - x 2 ki ■ k 3 
(1 - x 2 )k 1 ■ k 2 + a;ifc 2 2 ] [V (1) (l, 2, 3) + ^ (1) (3, 2, 1)1 tp {2 \l, 2, 3) 
+ (fei x fe 2 )[^(l,2,3)(fe 3 x V 3 )^ (2) (l,2,3) 
-^(1,2,3)^x^)^(1,2,3) 
- 3, 2) (fe 3 x V s ) (V 2) (l, 3, 2) - ^ (2) (2, 3, 1)) 

+ V (2) (l, 3, 2) (fc s x V 8 ) 3 ' 2 ) + <A (1) (2, 3, 1))] } . 

In Eq. (140 p we used the following definition 

V^«(l, 2, 3) = [(1 - arj) V fel - x 2 V fc2 ] 2, fc 3 = -fci - fc 2 ) , (42) 

and similarly for V 2 . Furthermore, in Eq. (1411) the operator V 3 is defined as 

V 3 V (i) (l, 2, 3) = - (a* V fcl + x 2 V fc2 ) 2, fc 3 = -k x - k 2 ) . (43) 

Using the momentum conservation constraint ki + k 2 + fc 3 = 0, one finds £^' lz=0 = —£f lz=0 . 
One then recovers the fact that the total contribution from the l z — component is equal 
to zero 

0. 



(44) 



For the l z = 1 component, we find: 
for the total u-quark contribution 

£ u z ' lz=1 = J [dx] 3 [d 2 fc] 3 {[fc 1 -fc 2 ^(l,2,3) + fc 2 V 4 )(l,2,3) 

- 2x^(1, 2, 3) + 2(1 - x 2 )^\l, 2, 3) - x^ 3 \l, 3, 2) - x 3 ^(l, 3, 2) 
feiV^Cl, 2, 3)) + fei • fc 2 ^ (4) (l, 2, 3)1 [(1 - ^^(l, 2, 3) - x 2 ^ (4) (l, 2, 3) 



x 



+ 



fc 1 -fc 3 ^ 3 )(l,2,3) + fc 2 -fc 3 ^( 4 )(l,2,3) 

x^ (3) (l, 2, 3) + x 2 ^ (4) (l, 2, 3) + x^ (3) (l, 3, 2) - (1 - x 3 )^\l, 3, 2) 

- (fei x fc 2 )^ (3) (l, 2, 3) [(fc 2 x V 2 + fc 3 x V 3 V 4) (1, 3, 2) - (fc 2 x V 2 )^(l, 2, 3) 

- (fc x x fc 2 )^ 4 )(l,2,3) [(fc 2 x V 2 + fc 3 x V 3 ) (^ (3) (1,3,2)-V {4) (1,3,2)) 
fc 2 xv 2 )^(l,2,3)l) ; (45) 



11 



for the total <i-quark contribution 



nd,l z =l 



[dz] 3 [d 2 k] 3 | \kfy® (1, 2, 3)) + fc x ■ fc 2 ^ 4 ) (1, 2, 3) 
(1 - X!)^ (3) (l, 2, 3) - x 2 ^ (4) (l, 2, 3) + (1 - X!)^ (3) (l, 3, 2) - x 3 ^ (4) (l, 3, 2) 
fc x ■ fc 3 ^ 3) (l, 2, 3) + fc 2 • k 3 ^ 4 \l, 2, 3)] [x^( 3 )(l, 2, 3) + x 2 ^ (4) (l, 2, 3) 
-(fci x fc 2 )^ (3) (l,2,3) [(fei x Vx)^ 4 )(l,3,2)- (kt x Vx + fca x V 3 )^ (4) (l, 2, 3) 
- (k, x fc 2 )^ (4) (l,2,3) [(fc! x V x ) (^(1,3,2) -^W(l,3,2)) 

+ (fei x Vi + fc 3 x V 3 )^ (3) (l, 2, 3)] } . (46) 
The sum of the u and d contributions from the l g — 1 component gives 
' [dx] 3 [d 2 fc] 3 

2^ (3) (1, 2, 3) + ^ {3) (1, 3, 2)1 \kfy®(l, 2, 3) + fci • fc 2 ^ (4) (l, 2, 3) 



+ 2^(1, 2, 3) [fei • fe 2 V (3) (l, 2, 3) - fc 2 V (4) (l, 2, 3) 
+ V (4) (l, 3, 2) [k, • M' (3) (l, 2, 3) + fc 2 ■ fc 3 ^ (4) (l, 2, 3) 
= Ph=i- 

For l z = —1 component, we find: 
for the total w-quark contribution 

C'' =_1 = / [dx] 3 [d 2 fc] 3 {[fc 1 -M' (5) (l,2,3)-fc 1 -fc 3 ^(l,3,2) 
x x 2 ^ (5) (l, 2, 3) - x 3 ^ 5 \l, 3, 2) - x 3 ^ {5) (2, 3, 1) - (1 - x 1 )^ 5 \2, 1, 3) 
+ [fe^ (5) (l,2,3)-fe 2 -fc 3 ^ (5) (l,3,2) 

x x 1 ^ (5) (2, 1,3)- x 3 ^ 5 \l, 3, 2) - x 3 ^ {5) (2, 3, 1) - (1 - x 2 )^ 5 \l, 2, 3) 
+ (fe 2 x fc 3 ) [^( 5 )(1, 3, 2) (fei x Vx + k 2 x V 2 ) (> 5) (2, 1, 3) - ^ 5 >(1, 2, 3) 
+ V (5) (l> 2, 3) ffej x Vi + k 2 x V 2 ) (V (5) (l, 3, 2) + ^(2, 1, 3) + ^ (5) (2, 3, 1) 



(47) 



(48) 
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for the total rf-quark contribution 

= J [dx] 3 [d 2 fc] 3 {[fc 2 • fc 3 ^ (5) (l,2,3) - fc 3 V 5 )(l,3,2) 

x ^ (5) (2, 1, 3) + x 2 ^ (5) (l, 2, 3) + (1 - x 3 ) (V 5) (l, 3, 2) + V (5) (2, 3, 1) 
+ (fe 2 x fc 3 ) ^ (5) (1, 3, 2) (k 3 x V 3 ) (V {5) (2, 1, 3) - ^ 5 >(1, 2, 3)) 
+ V (5) (l, 2, 3) (fc 3 x V 3 ) fv (5) (l> 3, 2) + ^ (5) (2, 1, 3) + ^ (5) (2, 3, 1) 



Adding the u and d contributions, we find that the l z 
total orbital angular momentum as 



(49) 



1 partial wave contributes to the 



"2 



+ 



[dx] 3 [d 2 fc] 3 [ - [fci • fc 2 ^(l,2,3) - fej ■ fc 3 ^( 5 )(l,3,2)]^ (5) (2,l,3) 

fc^ (5) (l, 2, 3) - fe 2 • M> {5) (1, 3, 2)] ^(1, 2, 3) 
fc 2 • M> (5) (1, 2, 3) - fc 3 2 V (5) (l, 3, 2)] [^ (5) (1, 3, 2) + ^ (5) (2, 3, 1) 
Pj,=-i- 



(50) 



For l z — 2 component, we find: 
for the total u-quark contribution 



nn,/ z =2 



[dx] 3 [d 2 fc] 3 4 fc 2 (1 - xi)kl - x 3 fci • fc s 



x 



?/; (6) (l, 2, 3) f2^ (6) (l, 2, 3) + ^ (6) (3, 2, 1)) + 7/^(3, 2, 1) (V 6) (l, 2, 3) + ^ (6) (3, 2, 1) 



fc 2 



x 2 fc. 2 + x 3 fc 2 • k c 



X 



x 



^ (6) (1, 2, 3) (> 6 )(2, 1,3)- ^ 6 >(3, 1, 2)) - ^ {6) (3, 2, 1)^ {6) (3, 1, 2) 



A; 2 (1 - xi)k 2 ■ k 3 - x 2 k x ■ k 3 



x 



^(1, 2, 3) (2^(1, 3, 2) + ^(2, 3, 1)J + ^ (6) (3, 2, 1)^ (6) (1, 3, 2) 
zjfei • fc 2 fc 2 + x 3 fc 2 • fc 3 fc 2 ] ^ 6 >(1, 2, 3)^ (6) (1, 2, 3) 
(1 - x 2 )fc 2 ■ fc 3 fc 2 - x 1 (2k 1 ■ k 2 k 2 ■ k 3 - kf fc 3 fc 2 2 ) 

^ {6) (1, 2, 3) (V 6 )(l, 3, 2) + ^ (6) (2, 3, 1)) - ^ 6 >(3, 2, 1)^ (6) (2, 3, 1) 



+ {ki x fc 2 )^ C6) (l,2, 3)[fc|(fe! x Vx + fc 2 x V 2 )^(2,l,3) 
-fc 3 2 (fe 1 xv 1 + fc3xV 3 )^(3 1 l,2) 

+fc 2 (fc 1 x V 1 )^(l,3,2) + fc 2 (fc 1 x Vi + fc 2 x V 2 )^(2,3,l)]} ; 



(51) 
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for the total <i-quark contribution 



0dJ,=2 



[dx] 3 [d 2 fc] ; 



xifei • k 2 kl + x 3 k 2 ■ fc 3 fc 2 



^ (6) (1,2,3) ^ (6) (1,2,3) + ^ (6) (3,2,1) 



[1 - x 2 )k 1 ■ k 2 kl - x 3 (2ki ■ k 2 k 2 ■ k 3 - k 1 ■ k 3 kr 



fc 2 



k 



x ^(3, 1, 2) ^(1, 2, 3) + ^ (6) (3, 2, 1) 

x 2 k, • fc 3 - (1 - z 3 )fci ■ fea] ^ (6) (1, 2, 3)^(2, 1, 3) 
Xlkl ■ k 3 - (1 - x 3 )fc 2 ] ^ 6 >(1, 2, 3)^ (6) (1, 2, 3) 
f xxfei • k 2 + x 2 kf\ ^ 6 \l, 2, 3) U (fl) (l, 3, 2) + ^ (6) (2, 3, 1) 



fei x fc 2 )^ (6) (l,2,3) [fcf(fc 2 x V 2 )^(3,l,2) + fc 2 (fe 3 x V 3 )^ {6) (2,1,3) 



fc 2 (fc 2 x V 2 + k 3 x V 3 )^ (6) (l, 3, 2) + fc 2 (fc 3 x V 3 )v> (6) (2, 3, 1) 



(52) 



The sum of the u and d contributions in the l z — 2 component to the total orbital angular 
momentum is 



X 



[d*] 3 [d 2 fc] 3 {^(1,2,3) 

2fc 2 fc 2 (V 6) (l, 2, 3) + ^ (6) (3, 2, 1)) - 2fc 2 fc 1 • fc 2 (V 6) (3, h 2) 
- 2fc 2 fc 2 ■ k 3 (2^ 6 >(1, 3, 2) + ^ (6) (2, 3, 1))] } 
2p/,= 2 . 



^ (6) (2,1,3) 



(53) 



V. RESULTS IN LIGHT-CONE QUARK MODELS 

The method developed in the previous sections is applied here within two light-cone quark 
models, a light-cone constituent quark model (LCCQM) 22|-|26[ and the light-cone version of 
the chiral quark-soliton model (LCxQSM) restricted to the three-quark sector 27H29|. These 
two models were recently applied to describe the valence-quark structure of the nucleon as 
observed in parton distribution functions, like generalized parton distributions, transverse- 
momentum dependent parton distributions, and form factors of the nucleon, giving a typical 
accuracy of about 30% in comparison with available data in the valence region. 
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TABLE I: The results for the quark orbital angular momentum from the LCCQM and the LCxQSM for 
u-, d- and total (u + d) quark contributions. 



Model 


LCCQM 




LCxQSM 




q 


u d 


Total 


u d 


Total 


it Eq. (SO]) 


0.131 -0.005 


0.126 


0.073 -0.004 


0.069 



In the LCCQM the nucleoli state is described by a LCWF in the basis of three free 
on-shell valence quarks. The three-quark state is however not on-shell Mjv 7^ .Mo = °°i 
where Ui is the energy of free quark % and is the physical mass of the nucleon bound 
state. The nucleon wave function is assumed to be a simple analytic function depending 
on three free parameters (including the quark mass) which are fitted to reproduce at best 
some experimental observables, like e.g. the anomalous magnetic moments of the proton 
and neutron and the axial charge. 

In the LCxQSM quarks are not free, but bound by a relativistic chiral mean field (semi- 
classical approximation). This chiral mean field creates a discrete level in the one-quark 
spectrum and distorts at the same time the Dirac sea. It has been shown that the distortion 



271 ] . Even though the 



can be represented by additional quark-antiquark pairs in the baryon 
xQSM naturally incorporates higher Fock states, we restrict the present study to the 3Q 
sector. 

Despite the apparent differences between the LCCQM and the xQSM, it turns out that 
the corresponding LCWFs have a very similar structure (for further details, we refer to [12 



301] ). The corresponding predictions from the LCCQM and the LCxQSM for u-, d- and total 
(u + d) quark contributions to the OAM are reported in Table [H We note that there is more 
net quark OAM in the LCCQM (J2 q ^l = 0.126) than in the LCxQSM = 0.069). For 

the individual quark contributions, both the LCCQM and the LCxQSM predict that l q z are 
positive for u quarks and negative for d quarks, with the w-quark contribution larger than 
the <i-quark contribution in absolute value. These results refer to the low hadronic scales 
of the models. Before making a meaningful comparison with results from experiments or 
lattice calculations which are usually obtained at higher scales, a proper treatment of the 
effects due to QCD evolution is essential 311 ] . However, this is beyond the scope of the 
present paper. 
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The explicit calculation within the LCCQM of the wave-function amplitudes in 
Eqs. (I36ap - fl36dj) can be found in Ref. 24|. Using these results and the expressions in 
the appendix, we can also calculate the distribution in x of the OAM £|, separating the 
contribution from each partial wave. The corresponding results for u and d quarks are 
shown in Fig. [TJ The x dependence of the different partial-wave amplitudes is very similar 




FIG. 1: Results for the distribution in x of the OAM £ z in the proton for u (left) and d (right) 
quark. The curves correspond to the contribution of the different partial waves: long-dashed curves 
for the light-cone amplitude with l z = 0, dotted curves for the light-cone amplitude with l z = 1, 
dashed-dotted curves for the light-cone amplitude with l z = —1, and short-dashed curves for light- 
cone amplitude with l z = 2. The solid curves show the total results, sum of all the partial-wave 
contributions. 

for u and d quarks. However, the total contribution has a distinctive behavior for u and 
d quarks, coming from a quite different interplay between the different partial-waves. For 
the u quarks, the dominant contribution comes from the l z = 1 amplitude (dotted curve), 
with positive sign, while the positive contributions coming from the l z — (long-dashed 
curve) and l z = 2 (short-dahsed curve) amplitudes are largely compensated by the negative 
contribution coming from the l z = — 1 amplitude (dashed-dotted curve). For the d quarks, 
the OAM arises from the competition between the positive l z = 1 and l z = 2 contributions, 
and the negative l z — — 1 and l z = contributions. As a result, the OAM for d quarks is 
much smaller than for u quarks, and goes from negative to positive values at x ~ 0.3. 
The integral over x of the different distributions in Fig[T]gives the value for u- and <i-quark 
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TABLE II: Results from the LCCQM for the contribution of the different partial waves to the OAM 1%. 
The first and second rows show the values for the u and d quarks, respectively, while in the third row are 
reported the sum of the u- and d-quark contributions. In the last row are the results for the squared norm of 
the different partial-waves, giving the probability to find the proton in a three-quark state with eigenvalue 
L of total OAM. 





lz = 


h = 1 


lz = -1 


l z = 2 


Total 


pu 


0.013 


0.139 


-0.046 


0.025 


0.131 


£i 


-0.013 


0.087 


-0.090 


0.011 


-0.005 


4 





0.226 


-0.136 


0.036 


0.126 


Pi, 


0.620 


0.226 


0.136 


0.018 


1 



OAM reported in Table [TTJ In the last row we also show the results for the squared norm 
of the different partial- waves pi z , giving the probability to find the proton in a three-quark 
state with eigenvalue l z of total OAM, according to Eq. (13"3"j) . 
It is interesting to rewrite the expression ffTSl) for the OAM as 

P M = J d 2 b(bx(k)i) z , (54) 

where (k) q is the distribution in impact-parameter space of the quark mean transverse 
momentum 

(fe)9(b) = Jdx d 2 k k p [ ll ]q (b, fe, x). (55) 

This distribution is shown in Fig. [2] for both u and d quarks. First of all, it appears that 
the mean transverse momentum (k) q is always orthogonal to the impact-parameter vector 
b. This is not surprising since a nonvanishing radial component of the mean transverse 
momentum would indicate that the proton size and/or shape are changing. The figure also 
clearly shows that u quarks tend to orbit anticlockwise inside the nucleon, corresponding to 
l u z > since the proton is represented with its spin pointing out of the figure. For the d 
quarks, we see two regions. In the central region of the nucleon, |b| < 0.3 fm, the d quarks 
tend to orbit anticlockwise like the u quarks. In the peripheral region, |6| > 0.3 fm, the 
d quarks tend to orbit clockwise. All this is consistent with the three-dimensional picture 
provided by the generalized parton distributions which indicates that the central region is 
dominated by the large x values, while the peripheral region is dominated by the low x 
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values. The approximate cancellation between the central (large x) and peripheral (small 
x) contributions leads then to a very small value for the <i-quark OAM. 



0.6 
0.4 
0.2 

's 

fa o.o 
-cT 

-0.2 
-0.4 
-0.6 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 




i 



FIG. 2: Distributions in impact parameter space of the mean transverse momentum of unpolarized 
quarks in a longitudinally polarized nucleon. The nucleon polarization is pointing out of the plane, 
while the arrows show the size and direction of the mean transverse momentum of the quarks. The 
left (right) panel shows the results for u (d) quarks. 



<jfc y > u [Gev/fni 2 ] 
0.2^ r\ 




bjfm] 



FIG. 3: Results for the different partial- wave contributions to the u quark mean transverse mo- 



mentum (k) u 



z y ) u e y as a function of b = b x e 2 



Figs. [3] and |4] show the different partial-wave contributions to the u and d quark mean 
transverse momentum (k) q as a function of b. Due to the axial symmetry of the system, 
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</t y > d [Gev/fm 2 ] 
0.3 r 



/ oi ; 



-1.5 



L=l 



L=0 



-1.0 




-0.3 



1.5 



b x [fm] 



L=- 



L=2 



Total 



FIG. 4: Results for the different partial- wave contributions to the d quark mean transverse mo- 
mentum (k) d = {k y ) d e y as a function of b = b x e x . 

it is sufficient to plot (k) q = (k y ) q e y as a function of b — b x e x . Similarly to Fig. [U for 
the u quarks, the dominant contribution comes from the l z = 1 amplitude (dotted curve), 
with positive sign, while the positive contributions coming from the l z = (long-dashed 
curve) and l z = 2 (short-dahsed curve) amplitudes are largely compensated by the negative 
contribution coming from the l z = —1 amplitude (dashed-dotted curve). For the d quarks, 
the OAM arises from the competition between the positive l z = 1 and l z = 2 contributions, 
and the negative l z — — 1 and l z = contributions, with a delicate balance between the 
different partial-wave contributions. As a result the total OAM takes small positive values 
at small \b\, becomes slightly negative at |b| « 0.3 fm, and vanishes at the periphery. 

VI. CONCLUSIONS 

In summary, we derived the relation between the quark OAM and the Wigner distribution 
for unpolarized quark in a longitudinally polarized nucleon. This relation is exact as long 
as we neglect the contribution of gauge fields, and provides an intuitive and simple repre- 
sentation of the quark OAM which resembles the classical formula given by the phase-space 
average of the orbital angular momentum weighted by the density operator. We compare 
this derivation with the LCWF representation of the OAM. The advantage in using LCWFs 
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is that they are eigenstates of the total OAM for each A-parton configuration in the nu- 
cleon Fock space. As a consequence, the total OAM can be simply calculated from the sum 
of squared LCWFs multiplied by the corresponding eigenvalues of the OAM operator. In 
the three-quark sector, we further decomposed the nucleon state in different partial-wave 
amplitudes, calculating the corresponding contributions to the quark OAM. 

These two representations of the OAM are equivalent, and allow to visualize comple- 
mentary aspects of the orbital motion of the quarks inside the nucleon. As examples, we 
adopted two different light-cone quark models, and discussed the corresponding results for 
the distribution in x of the OAM, as obtained from the LCWF overlap representation, as 
well as the distribution of the mean transverse momentum in the impact-parameter space, 
as obtained from the Wigner distributions. 
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Appendix A: Partial-wave decomposition of the distributions in x of the OAM 

In this appendix we summarize the results for the partial-wave contributions to the distri- 
bution in x of the OAM. In particular, we separately list the results for £ q z ulz , corresponding 
to the contribution of the ith quark with flavor q in the Fock state component with total 
OAM l z . 

For the l z = component, we find: 
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for the u quark 



nui,l z =0. 



= J [dx] 3 [d 2 k] 3 5(x-x 1 ){- ^ (1) (1,2,3) 
2((1 - x 1 )k 1 ■ k 2 + x 2 k(\^ 2 \\, 2, 3) + (x 2 fcx ■ k 3 - xsfci • fe 2 )</> (2) (3, 2, 1) 
-(fei x fe 2 ) [^(1,2,3)^! x V 1 )(-2^ 2 )(l,2,3)+^ (2) (3,2,l)) 
+ ^ (2) (1, 2, 3) x Vi) (2^(1, 2, 3) + ^ (1) (3, 2, 1))] } , (Al) 



0U2,lz=0, 



= j [dz] 3 [d 2 fc] 3 5(x - a; 2 ){^ (1) (l, 2, 3) [(1 - x 2 )k l ■ k 2 + xife 2 ] ^ 2 \l, 2, 3) 

+ (fei x fc 2 ) 2, 3) (fe 2 x V 2 )^ (2) (l, 2, 3) - ^ {2) (1, 2, 3) (fc 2 x V 2 )^ (1) (l, 2, 3)] } , 

(A2) 

t z ^=\x) = j [dx] 3 [d 2 k] 3 5(x-x 3 ){- ^ (1) (1,2,3) 

x (safe! • fc 3 - xifc 2 ■ fe 3 )^ (2) (l, 2, 3) + ((1 - x 3 )fc 2 • fe 3 + x 2 fe 2 )^ (2) (3, 2, 1) 

- (fci x fc 2 ) 2, 3) (fc 3 x V 3 ) ( - ^ {2) (1, 2, 3) + ^ (2) (3, 2, 1)) 

+ ^ (2) (1, 2, 3) (fe 3 x V s ) 2, 3) + ^ (1) (3, 2, 1))] }; (A3) 



for the d quark 



gd,2,lz=0 



= J [dx] 3 [d 2 fe] 3 5(x-a: 2 ){^ (1) (l,2,3) 
(1 - x 2 )k 1 ■ k 2 + xifc 2 2 ) ^ 2 \l, 2, 3) + (a; 3 fc 2 + (1 - a; 2 )fe 2 • fe 3 )^ (2) (3, 2, 1) 
- (fci x fe 2 ) [^(l, 2, 3) (fe 2 x V 2 ) ( - ^ (2) (1, 2, 3) + ^ 2 >(3, 2, 1)) 
+ ^ 2 \l, 2, 3) (fe 2 x V 2 ) 2, 3) + ^ (1) (3, 2, 1))] }, (A4) 



,i*=o, 



= J [dx] 3 [d 2 k] 3 5(x-x 3 ){-(x 2 k 1 -k 3 -x 1 k 2 -k 3 ^ 1 \l,2,3)^ 2 \l,2,3) 

+ (fci x fc 2 ) 2, 3) (fc 3 x V 3 )^ (2) (l, 2, 3) - ^ {2) (1, 2, 3) (fe 3 x V 3 )^ (1) (l, 2, 3)] } . 

(A5) 
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For the l z = 1 component, we find: 
for the u quark 



(!U\,l z = l 



[dx] 3 [d 2 fc] 3 S(x — x\) 

kfy®{l, 2, 3)) + fe x • M (4) (l, 2, 3)] [(1 - x0^ (3) (l, 2, 3) - x 2 ^ (4) (l, 2, 3) 
(fe! x ^[^(l^)^ x Vi)^(l,2,3) - ^(l,^)^ x ^^(1,2,3)]} 



(A6) 



"?*'=\x) = J [dx] 3 [d 2 k] 3 S(x - x 2 )| [fex • fe 2 ^ (3) (l, 2, 3)) + fc 2 V 4) (l, 2, 3) 
x - 2x1^(1, 2, 3) + 2(1 - x 2 )^ 4 \l, 2, 3) - x^ (3) (l, 3, 2) - x 3 ^ (4) (l, 3, 2) 
- (fe x x k 2 ) [V 3) (l, 2, 3) (fc 2 xV 3 )(- 2^W(1, 2, 3) + 3, 2) 

+■ ^ (4) (1, 2, 3) (fc 2 x V 2 ) (2^(1, 2, 3) + ^ (3) (1, 3, 2) 
-^(l,2,3)(fe 2 xV 2 )^(l,3,2) 



(A7) 



(W3,l z =l. 



[dx} 3 [d 2 k] 3 5(x - x 3 ){ [fei • fc 3 ^ (3) (l,2,3)) + k 2 • M> (4) (M,3) 

- x^ (3) (l, 2, 3) - x 2 ^ (4) (l, 2, 3) - x 1 ^ 3 \l, 3, 2) + (1 - x 3 )^ (4) (l, 3, 2) 
- (fei x fe a ) [^ (3) (1, 2, 3) (fc 3 x V 3 ) ( - ^ (4) (1, 2, 3) + V (4) (i, 3, 2)) 
+ ^ (4) (1, 2, 3) (fcs x V s ) (> 3) (1, 2, 3) + ^(1, 3, 2)) 

-#>(l,2,3)(fc 3 x V 3 )#>(1,3,2)]} ; (A8) 



for the d quark 



,lz=l, 



[dx} 3 [d 2 k] 3 6{x - xOj [jfejty (3) (l, 2, 3)) + fei • fc 2 ^ (4) (l, 2, 3) 
(1 - X!)^ (3) (l, 2, 3) - x 2 ^ (4) (l, 2, 3) + (1 - X!)^ (3) (l, 3, 2) - x 3 ^ (4) (l, 3, 2) 
- (fci x k 2 ) [^ (3) (1, 2, 3) x Vi) ( - ^ (4) (1, 2, 3) + ^(1, 3, 2)) 
+ ^(1, 2, 3) (kt x Vx) (^ {3) (1, 2, 3) + ^ (3) (1, 3, 2)) 

-^(1,2,3)^x^)^(1,3,2))]} , (A9) 
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{ 



[dx] 3 [d 2 k] 3 8(x - x 3 ) 

fcx • fc 3 ^ (3) (l, 2, 3)) + k 2 ■ M> (4) (1, 2, 3)] [ - x^ (3) (l, 2, 3) - x 2 ^ (4) (l, 2, 3) 
(ki x fc 2 )[^ 3 )(l,2,3)(fc 3 x V 3 )^ (4) (l,2,3)) -V (4 )(l,2,3)(fc 3 x V 3 )^ (3) (l, 2, 3)] } 



For l z — — 1 component, we find: 
for the it quark 

x 2x 2 ^ (5) (l, 2, 3) - (1 - x0^ (5) (2, 1, 3) + x 2 ^ (5) (3, 2, 1) - 2x 3 ^ (5) (l, 3, 2) 
+ (1 - X!)^ (5) (3, 1,2)- x 3 ?/; (5) (2, 3, 1)1 - (fej x fc 2 )^ (5) (l, 2, 3) 



x ( fci x Vi 



'^ (5) (3, 1, 2) - ^ (5) (2, 1, 3) - 2^ 5 >(1, 3, 2) + ^(2, 3, 1) 



QU2,lz=— 1 



flU3,Z z = -l, 



[dx] 3 [d 2 fc] 3 5(x - x 3 ) I fc 2 • k^ (1,2,3) 

x 2 ^ (5) (l, 2, 3) + x 2 ^ 5 )(3, 2, 1) + (1 - x 3 )^ (5) (l, 3, 2) - x^ (5) (3, 1, 2) 
- (fci x fc 2 )^ 5 >(l, 2, 3) (fc 3 x V 3 ) U (5) (3, 1, 2)) - ^ 5 >(1, 3, 2) 



(A10) 



(All) 



[dx] 3 [d 2 k] 3 5(x - z 2 )|fc 2 fy( 5 >(l, 2, 3) 
(1 - x 2 )^ (5) (l, 2, 3) + x^ (5) (2, 1, 3) - x 3 ^ (5) (l, 3, 2) - x# (5) (2, 3, 1) 
(fci x fc 2 )^ 5) (l, 2, 3) (fc 2 x V 2 ) [^ (5) (2, 1, 3) + ^(1, 3, 2) + ^(2, 3, 1)] } , (A12) 



(A13) 



for the d quark 

g<h,i.=-i(x) = j [dx} 3 [d 2 k] 3 5(x-x 2 ){ -Jb 2 V (5) (l,2,3) 

x (1 - x 2 )V (5) (l, 2, 3) + (1 - x 2 )^ 5) (3, 2, 1) + x 3 ^ {5) (l, 3, 2) + x^ (5) (3, 1, 2) 
- (fci x fc 2 )V (5) (l, 2, 3) ffe a x V 2 ) U (5) (3, 1,2)- ^ (5) (1, 3, 2) 



(A14) 



e*» l ' = -\x) = j [dx] 3 [d 2 fc] 3 5(x-x 3 )|fc 2 -fc 3 ^ (5) (l,2,3) 

x Xl if>®{2, 1, 3) + x 2 if>®(l, 2, 3) + (1 - x 3 )^ (5) (l, 3, 2) + (1 - x 3 )^ {5) (2, 3, 1) 
+ (fc x x fe 2 )^ (6) (l, 2, 3) ffc 3 x V 3 ) U 5 >(2, 1, 3) + ^(1, 3, 2) + ^ 5 >(2, 3, 1) 



(A15) 



23 



For l z — 2 component, we find: 
for the u quark 



QUI, 



X 



+ 



f [dx] 3 [d 2 k] 3 S{x - Xl )^k^ 6 \l,2,3) 
(1 - Xl )kl - x 3 k x • fc 3 ] [2^ (6) (1, 2, 3) + ^ (6) (3, 2, 1) 



^(3,1,2) -^(2,1,3) 
(1 - Xl )k 2 ■ k 3 - x 2 k x • fe 3 ] [2^( 6 )(1, 3, 2) + ^ (6) (2, 3, 1) 
+ (fcx x fe 2 )^ (6) (l,2, 3)[fe|(fe! x Vi) ^ (6) (2,1,3))-^ (6) (3,1,2) 
+ k 2 (k, x V x ) (2^ (6) (1, 3, 2) + ^(2, 3, 1))] } , 



■I' 



[dx] 3 [d 2 k] 3 S(x - x 2 ) 

x x klki ■ k 2 + x 3 k 2 k 2 ■ k 3 ] ^ (6) (1, 2, 3) ^ (6) (1, 2, 3) 
(1 - x 2 )k 2 k 2 ■ k 3 — xi (2k x ■ k 2 k 2 ■ k 3 - k 2 k x ■ fc 3 j 

2, 3)^(1, 3, 2) + ^(2, 3, 1) 2, 3) - ^(3, 2, 1) 



(k.xk,)^ (1,2,3) 



k 3 (k 2 x V 2 



1,3) 



+ fc 2 (fc 2 x V 2 ) (^(l, 3, 2) + V> (6) (2, 3, 1))] } , 



X 



,lz=2. 



X) 



j [dx] 3 [d 2 k} 3 6(x - x 3 ){^ 6 \l,2,3)k 



(1 - x 3 )fc 2 - zifci • fe 3 ) (V 6) (3, 2, 1) + ^ (6) (1, 2, 3)) + (xifca ■ fc 2 + x 2 k 2 

- ((1 - x 3 )fci • k 2 - x 2 k x ■ fc 3 )^ (6) (3, 1, 2) 

- (fei x fc 2 )^ 6 >(l,2,3) [k 2 (k 3 x V 3 )^ (6) (3,l,2) -kl(k 3 x V 3 )^ {6) (1, 3, 2) 
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for the d quark 



,Z*=2/ 



\x 



H 3 [d 2 k]j( 



x - x 2 



x 1 k*k 1 ■ k 2 + x 3 fc 2 fc 2 • fe 3 l ^ (1, 2, 3) U {6) (1, 2, 3) + ^ (6) (3, 2, 1 
1 - x 2 )kfk 2 ■ k 3 - xi(2fci • fc 2 fc 2 ■ fc 3 - fci • fc 3 fc 2 5 
x V (6) (l, 3, 2) U 6 )(l, 2, 3) + ^ (6) (3, 2, 1) 
+ x fe 2 )V> (6) (l,2,3) feffca x V 2 )^ (6) (l,3,2)-fc 3 2 (fc 2 x V 2 W 6) (3, 1, 2) 



, (A19) 



X 



[dx] 3 [d 2 fc] 3 5(x-x 3 ){fc 3 V (6) (l,2,3) 
((1 - x 3 )Aq ■ fc 2 - x 2 fci ■ fc 3 )^ (6) (2, 1, 3) - (sxfci ■ fe 3 - (1 - x 3 )fei 2 )^ (6) (l, 2, 3) 



+ fxifci ■ fc 2 + x 2 fc 2 ) (V 6 >(1, 3, 2) + ^ (6) (2, 3, + (fei x fc 3 )^ (6) (l, 2, 3) 



x 



fc! 2 (fc 3 x V 3 )(V' {6) (1,3,2) + ^ (6) (2,3,1)) +fc 3 2 (fc 3 x V 3 )V' {6) (2,1,3)]} . (A20) 
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